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ABSTRACT 

In recent years, an increasing number of publications have been addressed to the 
pecuhar and mysterious pre-main sequence star V582 Mon, also known as KH 15D. 
This extraordinary T Tauri star, located in the young star cluster NGC 2264, appears 
as to be an eclipsing variable. In the present paper, we report a unique and self- 
consistent set of light curves in the blue and near-infrared bands, spanning a 15-year 
interval (epoch 1955-1970). Our photometric data show clearly the beginning of the 
eclipse stage occurred in early 1958 in the blue, and perhaps around four years later 
in the infrared. The light curve period turns out to be the same reported by recent 
observations (about 48.3 days), so that no evidence for a period change results. On the 
other hand, in our data the light curve shape appears as sinusoidal and is therefore 
different from the one displayed today. The photometric behaviour, determined with 
time-series and colour-index analysis, suggests that V582 Mon (KH 15D) could be 
initially surrounded by an accretion disk/torus seen edge-on, with subsequent thin 
dust formation at the beginning of the blue radiation absorption. The dust could then 
aggregate into larger particles providing the transition between selective and total 
absorption, accompanied with eclipsing variability in the infrared. The minima of the 
periodic light curve become deeper due to the increasing dimension and number of 
dust grains, and then flattens due to a contraction in the disk. 

Key words: stars: individual (V582 Mon) — stars: individual (KH 15D) — stars: 
variables: other — stars: pre-main sequence — stars: circumstellar matter — tech- 
niques: photometric — stars: planetary systems: protoplanetary disks 



1 INTRODUCTION 

The variable star V582 Mon {a = 6h 41m 10.3s, 5 = 
+09° 28' 34" (J2000.0)) was discovered by 
iBadalian fc Erastoval lll972f) . It is classified as a type 
INS (i.e. Orion variable with ra pid light chanRe s: up to 1 
mag in 1-10 days) in the GOVS iKholoDoylll98Rl^. This star 
was re -discovered as an eclipse variable bv lKearns 8z HerbstI 
l)l998h during a systematic study of young star variability 
in the region of the cluster NGC 2264. The s tar, named 
15D in a list of 209 objects teearns et al.lll997^ ■ appeared 
as an eclipsing variable with an amplitude of about 3 mag 
in the 7c (Cousins) band, with a period of 48-49 days. This 
was a unique cas e among the very y oung stars observed 
as discussed by iKearns et alJ ( |l222l)- The star is also 
designated as n. 15 in the list oflp Sk _gt al.l (|200n^ and as 
n. 391 in the list of lFlaccomio et al .1 ^1999^ ■ 
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Table 1. The number of Blue and IR photographic plates in the 
field of V582 Mon (KH 15D) from 1955 to 1970. 



Telescope 


N Blue 


Years 


N IR 


Years 


Reflector Loiano 0.60m 


70 


1957-59 


9 


1957-59 


Reflector Asiago 1.22m 


14 


1955-63 


5 


1956-63 


Schmidt Asiago 0.50/0. 40m 


52 


1959-62 


38 


1959-65 


Schmidt Asiago 0.90/0. 65m 


42 


1967-70 


47 


1967-70 


Schmidt Bergedorf 1.20/0. 80m 


5 


1962-64 







The depth (3.5 mag) and the long duration of the eclipse 
(about 40% of each cycle, which lasts for 48. 4 days) rule out 
a companion star as the cause of this event l|Hamilton et alJ 
l200ll:lHerbst et al.ll2002^ . Polarization levels of 2% observed 
at half eclipse an d 0.2 ± 0.2 observed out of the eclipse 
IjAgol et alJl2004l) . would imply that the obscuring mate- 
rial is composed of rather large particles. In order to explain 
these phenomena, two hypotheses were formulated: 1) a dou- 
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Figure 1. The complete light curves in the blue {~ B band) and photographic IR {~ Ic band) of V582 Mon (KH 15D) from 1955 to 
1970 (due to the largo plot scale nearby observations appear as a single point). 



ble star obscured in different ways by circumstellar matter; 
2) a single star surrounded by a peculiarly warped or ridged 
proto-planetary disk. In the second case, the light of the 
star is periodically intercepted by this dusty ridge orbiting 
around the star in a plane placed edge-on with respect to 
the line of sight. 

Actually, there is a large consensus that the eclipses 
are caused by circumstellar material, but it is not clear 
what is the composition and the spatial distribution of this 
material. Theories include an edge-on circumstellar disk 



with a mass den 
by a proto-plane1 


sity wave or a disk deformation induced 


iHamilton et al.ll200ll: 


iHerbst et al.ll2002t 


IWinn et alj 12003: Aeol et al. 




, an 


orbiting vortex of 



solid particles ijBarec & Viton 2003 ), and an asymmetric 
common envelope (Grinin fc Tambovtseval2002h . V582 Mon 
(KH 15D) is quite unique and extremely interesting be- 
cause the fortuitous ahgnment may allow us to study the 
circumstellar (or even protoplanetary) processes, and be- 
cause the eclipse provides a natural "coronagraph" that 
can be used to m ap out the environme n t of the underly- 
ing T Tauri star Hamilton et aP \2mA lAeol et al.l E)ol 
iDeming. Charbonneau. fc HarringtonI |2004|) . In fact, this 
object should be a weak-line T Tauri star, still accreting 
matter from a circums tellar disk and endowe d with a well- 
coUimated bipolar jet jHamilton et al.ll200j) . High resolu- 
tion near-IR spectroscopy revealed the presence of promi- 
nent molecular hydrogen emission probably due to shock 
excitation of the amb ient gas by a bipolar outflow from the 
circu mstellar disk (|Deming. Charbonneau. fc HarringtonI 
|2Qq3)- As ev idence of this hypothes is, there is also the ob- 
servation of iTokunaga et alJ i2004h , who revealed clearly 
an H2 emission fliament next to the central star. Fur- 
thermore, spectro-polarimetric observations suggest that 
the star is completely eclipsed, so that the flux during 



the eclipse is entirely due to scattered light. The scatter- 
ing should be due to large dust grains (~ 10 fim) sim- 
ilar to t he interplanetary grains that scatter the zodia- 
cal light jAgol et alj|2004lfl . Recent high-resolution spectro- 
scopic monitoring s howed that V582 M on is a single-line 
spectroscopic binary i Johnson et all2004l) . A m odel of grad- 
ual op aque-screen occultation was proposed bv lWinn et alJ 
l|2004) . and a model with an eccentric binary that is encir- 
cled by a vertically thin, inclined ring of dust and gas was 
suggested by (Chiang & Murray-Clav 2004) in order to ex- 
plain the available historical light curve prior to the present 
work. 

The light curve of V582 Mon (KH 15D) from 1967 to 
1982, based on the Asiago Observatory infrared plates, dis- 
plays a periodic alternation between brig ht and faint states 
with a constant period of 48.35 days djohnson fc Winnl 
l2004h . An international observation campaign during 2001- 
2002 showed clearly a periodic variabi lity recurring with a 
period of 48.36 days iHerbst e t al. 2002). Observations dur- 
ing six seasons from 1995 to 2001 showed again a periodical 
varia bility like an eclipsing binary, with a period of 48.34 
days iHamilton et alJl200J) . On the other hand, the analy- 
sis of the Harvard College Observatory archival plates from 
1913 to 1951 showed that there was no variability, and the 
deep eclipses obse rved today did not occur in the first half of 
the 20th century ijWinn et alJl2003t l. In that work, the need 
for further, more recent archival images was underlined, as 
they might have shown the onset of eclipses. 

In our paper, we present the results of the analysis of 
the plate images obtained by one of us (PM) during 15 years 
in a systematic study of the variability of young stars in the 
cluster NGC 2264. The data presented and discussed here 
pertain to the star V582 Mon (KH 15D) and concern the 
observations performed in the blue (B) and near-infrared 
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Table 2. Our photometric measurements in the blue and IR of V582 Mon (KH 15D) in the period 1955-1970 (for the meaning of the 
telescope acronyms see Table [Tl . 



l.lD-2430000] 



DiltL- 

l.lD-2430000] 















1 QflS (179'; 

lyoo.u i iO 






5427 549 


1955 8769 


16 3 


1 


RA 122 


6243 292 


1958 1103 


16 2 


1 


5507 438 


1956 0957 


16 3 


1 


RA 122 


6248 274 


1958.1240 


16 4 


1 


5508 384 


1956 0982 


16 3 


1 




6250 474 


1958 1300 


16 4 


1 








1 


RA 122 




1 Qf;S l/IRR 

iyoo.i4too 


16 4 


1 


5727 637 


1956 6985 


16 3 


1 


RA 122 


6275 367 


1958.1981 


16 4 


1 


5755 663 


1956 7753 


16 1 


1 




6277 441 


1958 2038 


16 4 


1 


5760.626 


1956 7S8S 


16 1 


1 


RA 122 


6285 347 


1958.2255 


16 3 


1 


5786 676 


1956 8602 


16 


1 




6303 336 


1958 2747 


16 6 


1 


5867 524 


1957 0815 


16 4 


1 


60^ 


6313 341 


1958 3021 


16 5 


1 


5868 531 


1957 0S43 


16 3 


1 


RA 122 


6314 326 


1958.3048 


16 3 


1 


5868 558 


1957 0843 


16 4 


1 




6330 340 


1958 3486 


15 9 


1 


5870 352 


1957 0893 


16 4 


1 


RL 60 


6460 607 


1958.7053 


16 2 


1 


5872 335 


1957 0947 


16 4 


1 


RL 60 


6579 406 


1959.0305 


16 1 


1 


5889 428 


1957 1415 


16 3 


1 




6595 240 


1959 0739 


16 3 


1 


5891 444 


1957 1470 


16 4 


1 


RL 60 


6596 253 


1959.0767 


16 5 


1 


5895 364 


1957 1577 


16 3 


1 




6598 292 


1959 0823 


16 5 


1 


5899 358 


1957 1687 


16 4 


1 


RL 60 


6599 434 


1959 0854 


16 5 


1 


5900.472 


1957 1717 


16 4 


1 


RL 60 


6600 352 


1959.0879 


16 3 


1 


5901 497 


1957 1745 


16 4 


1 




6606 350 


1959 1043 


16 6 


1 


5902 385 


1957 1770 


16 3 


1 


RL 60 


6607 384 


1959.1071 


16 4 


1 


OyU0.4:4:U 


1 Qt;7 1 7Qfi 




1 




OOUo.OZO 


lyoy.iuy / 


16 6 


1 


5918 336 


1957 2206 


16 5 


1 


RL 60 


6610 303 


1959 1151 


16 7 


1 


5946 348 


1957 2973 


16 4 


1 


RL 60 


6613 420 


1959.1237 


16 3 


1 


6083 605 


1957 6731 


16 5 


1 


RL 60 


6627 313 


1959 1617 


16 3 


1 








1 










1 


6104.631 


1957 7307 


16 '' 


1 


RL 60 


6630 381 


1959.1701 


16 5 


1 


6109 622 


1957 7443 


16 4 


1 




6636 343 


1959 1864 


16 3 


1 


6110 595 


1957 7470 


16 2 


1 


RL 60 


6643 330 


1959.2056 


16 1 


1 




1 Qt;7 7Q8Q 


16 2 


1 






lyoy.iuou 


16 3 


2 


6130.516 


1957.8015 


16.4 


0.1 


RL 60 


6661.317 


1959.2548 


16.2 


0.2 


6137.503 


1957.8207 


16.5 


0.1 


RL 60 


6662.292 


1959.2575 


16.2 


0.2 


6138.511 


1957.8234 


16.4 


0.1 


RL 60 


6670.323 


1959.2795 


16.5 


0.1 


6167.455 


1957.9027 


16.4 


0.1 


RL 60 


6812.615 


1959.6690 


16.0 


0.1 


6169.599 


1957.9086 


16.4 


0.1 


RL 60 


6818.610 


1959.6854 


15.9 


0.1 


6170.615 


1957.9113 


16.3 


0.1 


RA 122 


6846.622 


1959.7621 


16.5 


0.2 


6174.608 


1957.9223 


16.5 


0.1 


RL 60 


6952.351 


1960.0516 


16.3 


0.2 


6191.469 


1957.9684 


16.4 


0.1 


RL 60 


7016.400 


1960.2270 


16.3 


0.2 


6194.495 


1957.9767 


16.2 


0.1 


RL 60 


7259.499 


1960.8925 


16.4 


0.2 


6195.392 


1957.9792 


16.3 


0.1 


RL 60 


7319.515 


1961.0569 


16.2 


0.2 


6198.468 


1957.9876 


16.4 


0.1 


RL 60 


7336.387 


1961.1030 


16.5 


0.2 


6212.279 


1958.0254 


16.7 


0.1 


RL 60 


7341.333 


1961.1166 


16.6 


0.2 


6214.287 


1958.0309 


16.5 


0.1 


RL 60 


7345.272 


1961.1274 


16.5 


0.2 


6219.353 


1958.0448 


16.5 


0.1 


RL 60 


7561.595 


1961.7196 


16.8 


0.2 


6225.316 


1958.0611 


16.4 


0.1 


RL 60 


7562.594 


1961.7224 


16.5 


0.2 


6228.369 


1958.0695 


16. 4 


0,1 


RL 60 


7563.626 


1961.7252 


16.8 


0.2 



[.lU-2430000] 



Date 
[.lD-2430000] 



■»g] [mag] 



7564.633 


1961 7279 






SA 50/40 


8111 333 


1963 2247 






SB 120/80 


7619 581 


1961 8784 


16 6 


2 


SA 50/40 


8432 375 


1964 1037 


163 


1 


SB 120/80 


7621.639 


1961.8840 


16.7 


0.2 


SA 50/40 


9765.626 


1967.7539 


15.9 


0.1 


SA 90/65 


7641.403 


1961.9381 


16.4 


0.2 


SA 50/40 


9769.626 


1967.7649 


16.0 


0.1 


SA 90/65 


7642 590 


1961 9414 






SA 50/40 


9771 597 


1967 7703 






SA 90/65 


7645.423 


1961 9491 


16 4 


2 


SA 50/40 


9772 542 


1967 7729 


15 8 


1 


SA 90/65 


7648.488 


1961.9575 


16.4 


0.2 


SA 50/40 


9773.567 


1967.7757 


15.8 


0.1 


SA 90/65 


7649.515 


1961.9603 


16.4 


0.2 


SA 50/40 


9774.650 


1967.7786 


16.0 


0.1 


SA 90/65 


7651.613 


1961.9661 


16.3 


0.2 


SA 50/40 


9775.606 


1967.7813 


15.9 


0.1 


SA 90/65 


7659 306 


1961 9871 






SA 50/40 


9796 593 


1967 8387 






SA 90/65 


7671.496 


1962.0205 


16.7 


0.2 


SA 50/40 


9826.689 


1967.9211 


16.3 


0.1 


SA 90/65 


7672.461 


1962.0232 


16.8 


0.2 


SA 50/40 


9827.677 


1967.9238 


15.9 


0.1 


SA 90/65 


7680.592 


1962.0454 


16.3 


0.2 


SA 50/40 


9830.653 


1967.9320 


16.3 


0.1 


SA 90/65 


7689.300 


1962 0693 






SA 50/40 


9831 672 


1967 9348 






SA 90/65 


7692.306 


1962.0775 


16.3 


0.2 


SA 50/40 


9852.603 


1967.9921 


17.0 


0.1 


SA 90/65 


7694.235 


1962.0828 


16.3 


0.2 


SA 50/40 


9860.556 


1968.0138 


16.3 


0.1 


SA 90/65 


7695.401 


1962.0860 


15.9 


0.2 


SA 50/40 


9864.550 


1968.0248 


16.2 


0.1 


SA 90/65 


7696 293 


1962 0884 






SA 50/40 


9875 351 


1968 0543 






SA 90/65 


7698.292 


1962.0939 


16.4 


0.2 


SA 50/40 


9876.306 


1968.0570 


15.8 


0.1 


SA 90/65 


7699.294 


1962.0966 


16.3 


0.2 


SA 50/40 


9877.311 


1968.0597 


15.9 


0.1 


SA 90/65 


7700.346 


1962.0995 


16.4 


0.2 


SA 50/40 


9878.497 


1968.0630 


15.9 


0.1 


SA 90/65 


7701.307 


1962 1021 






SA 50/40 


9879 322 


1968 0652 






SA 90/65 


7702 292 


1962 1048 


16 7 


2 


SA 50/40 


9881 492 


1968 0712 


163 


1 


SA 90/65 


7717.274 


1962.1459 


16.8 


0.2 


SA 50/40 


9885.510 


1968.0822 


16.5 


0.1 


SA 90/65 


7719.306 


1962.1514 


16.4 


0.2 


SA 50/40 


9887.488 


1968.0876 


16.5 


0.1 


SA 90/65 


7720 310 


1962 1542 






SA 50/40 


9905 422 


1968 1367 






SA 90/65 


7731.321 


1962 1843 


15 9 


2 


SA 50/40 


9941 353 


1968 2351 


16 8 


1 


SA 90/65 


7732.285 


1962.1870 


16.3 


0.2 


SA 50/40 


9942.378 


1968.2379 


16.8 


0.1 


SA 90/65 


7752.316 


1962.2418 


16.2 


0.2 


SA 50/40 


9943.319 


1968.2404 


16.8 


0.1 


SA 90/65 


7755.309 


1962 2500 






SA 50/40 


9944 317 


1968 2432 






SA 90/65 


7757.310 


1962.2555 


16.4 


0.2 


SA 50/40 


9947.325 


1968.2514 


16.8 


0.1 


SA 90/65 


7760.356 


1962.2638 


16.5 


0.2 


SA 50/40 


10241.498 


1969.0568 


16.4 


0.1 


SA 90/65 


7762.322 


1962.2692 


16.4 


0.2 


SA 50/40 


10264.266 


1969.1191 


16.1 


0.1 


SA 90/65 


7939.597 


1962.7545 


16.1 


0.1 


RA 122 


10273.295 


1969.1439 


16.7 


0.1 


SA 90/65 


7972.562 


1962.8448 


16.2 


0.1 


SB 120/80 


10274.456 


1969.1470 


16.7 


0.1 


SA 90/65 


7974.564 


1962.8503 


16.3 


0.1 


SB 120/80 


10504.663 


1969.7773 


16.4 


0.1 


SA 90/65 


7991.583 


1962.8969 


16.3 


0.2 


SA 50/40 


10510.643 


1969.7937 


16.3 


0.1 


SA 90/65 


7993.556 


1962.9023 


16.4 


0.2 


SA 50/40 


10566.520 


1969.9467 


16.9 


0.1 


SA 90/65 


7998.497 


1962.9158 


16.0 


0.2 


SA 50/40 


10566.571 


1969.9468 


16.9 


0.1 


SA 90/65 


8001.487 


1962.9240 


16.5 


0.2 


SA 50/40 


10569.530 


1969.9549 


17.1 


0.1 


SA 90/65 


8003.556 


1962.9297 


16.8 


0.2 


SA 50/40 


10588.493 


1970.0068 


16.0 


0.1 


SA 90/65 


8020.522 


1962.9761 


16.4 


0.2 


SA 50/40 


10589.562 


1970.0098 


16.0 


0.1 


SA 90/65 


8046.310 


1963.0467 


16.5 


0.1 


RA 122 


10621.527 


1970.0973 


16.5 


0.1 


SA 90/65 


8049.480 


1963.0554 


16.4 


0.1 


SB 120/80 


10626.344 


1970.1105 


16.8 


0.1 


SA 90/65 


8050.266 


1963.0575 


16.3 


0.1 


RA 122 


10647.369 


1970.1680 


16.0 


0.1 


SA 90/65 


8057.410 


1963.0771 


16.4 


0.1 


RA 122 













EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
EL 60 
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SA 50/40 
SA 50/40 
SA 50/40 
EL 60 
EL 60 
EL 60 
SA 50/40 
SA 50/40 
SA 50/40 
SA 50/40 
SA 50/40 
SA 50/40 
SA 50/40 
SA 50/40 
SA 50/40 
SA 50/40 
SA 50/40 



Date 
[JD-2430000] 


Year 


I 

[mag] 


AI 

[mag] 


Tel. 


Date 
[.ID-2430n0fl] 


Year 


I 

[mag] 


AI 

[mag] 


Tel. 


Date 
[JD-2430000] 


Year 


I 

[mag] 


AI 

[mag] 


Tel. 


Date 
[,ID-2430000] 


Year 


I 

[mag] 


AI 

[mag] 


Tel, 


5755.645 


1956.7752 


13.9 


0.1 


RA 122 


7364.304 


1961-1795 


13.9 


0.2 


SA 50/40 


8050.281 


1963.0576 


14.6 


0.1 


BA 122 


9881.506 


1968,0712 


13,9 


0,1 


SA 90/65 


5815.572 


1956.9393 


14.0 


0.1 


BA 122 


7,365.352 


1961.1823 


13.9 


0.2 


SA 50/40 


8057.392 


1963.0770 


14.6 


0.1 


BA 122 


9885.525 


1968.0822 


14,3 


0.1 


SA 90/65 


6101.603 


1957.7224 


13.9 


0.1 


BL 60 


7366.372 


1961.1851 


13.8 


0.2 


SA 50/40 


8799.463 


1965.1087 


14.0 


0.2 


SA 50/40 


9887.503 


1968.0876 


14,4 


0.1 


SA 90/65 


6109.602 


1957.7443 


13.9 


0.1 


BL 60 


7,367.306 


1961.1877 


13.9 


0.2 


SA 50/40 


9763.640 


1967.7485 


13.9 


0.1 


SA 90/65 


9905.462 


1968.1368 


14,0 


0,1 


SA 90/65 


6130.538 


1957.8016 


14.0 


0.1 


BL 60 


7368.326 


1961.1905 


13.8 


0.2 


SA 50/40 


9765.647 


1967.7540 


14.1 


0.1 


SA 90/65 


9941.338 


1968.2350 


14,8 


0,1 


SA 90/65 


6138.535 


1957.8235 


13.9 


0.1 


BL 60 


7,369.30!) 


1961.1932 


14.0 


0.2 


SA 50/40 


9769.651 


1967.7650 


13.9 


0.1 


SA 90/65 


9942.358 


1968.2378 


14,8 


0,1 


SA 90/65 


6170.507 


1957.9110 


14.1 


0.1 
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(7) bands between October 13, 1955 and March 1, 1970. 
With these data, it is possible to determine when and how 
the eclipse phenomenon occurred, its first evolution, and the 
subsequent dramatic variations. 

In Section]^ we describe the telescopes used, the photo- 
graphic plates obtained, and the data reduction. In Section 

01 we show results from the stationary period up to 1970, 
and in Section |1] we present a brief discussion and our con- 
clusions. 

2 OBSERVATIONS 

The field of the young cluster NGC 2264 was photographed 
in the winter season 1955-56 with the 0.6m Loiano (Bologna, 



Italy) and the 1.2m Asiago (Vicenza, Italy) telescopes. 
These observations led to the discovery of 33 new vari- 
able stars, studied jointly with the other known 78 vari- 
able stars (Rosino. Grubissich, & Maffci 1957). In this sys- 
tematic study, V582 Mon (KH 15D) was not detected as a 
variable star. 

In the subsequent years, the search was extended by 
PM to the neighboring region of this cluster. In 15 years, 
183 blue and 99 photographic infrared plates were obtained, 
with five different telescopes. TableQshows the distribution 
of the photographic material collected. The number of plates 
used does not include some bad plates (e.g. scratched, or 
partially lighted, or with a covering halo due to the bright 
star in the field), and it does not comprise few underexposed 
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plates (in which the magnitude of the variable star was much 
fainter than the limit-magnitude of the plate). Hence these 
latter plates (about ten in B) were not useful to set meaning 
upper limits. 

The merit of the material used in this work lays in its 
uniformity. In the blue, an emulsion 103a-O without filter 
was used. In the photographic-infrared, an emulsion I-N hy- 
persensitized with an RG5 filter was used. The I-N/RG 5 
combination transmits the spectral range 6800-8800 A, and 
has a transmission curve similar to the 7000-9000 A band- 
pass of the Cousins Ic band. Moreover, a large fraction of 
the plates (about 60%, see Tab.^and Tab.|5J were obtained 
with the 0.60m Loiano reflector and the 0.90/0. 65m Asiago 
Schmidt telescope, having the same focal length (210 cm) 
and almost the same focal ratio (f/3.50 and f/3.23 respec- 
tively). 

The irregular variability of the studied stars required a 
high frequency in the plate collection. As a consequence, for 
V582 Mon (KH 15D), a series of long duration observations 
(except for the seasonal, weather and full-moon interrup- 
tions) are available, in the years when the star began to 
show its variability and during the first evolution phases. 
For the first years this is true especially for the blue plates; 
for the subsequent years also in the infrared ones. 

The magnitudes of comparison stars were tak en from 
the bl ue photoelectric and photographic sequences of lWalkeH 
il956h . after removing stars that appeared variable after- 
wards. The calibrations curves were constructed using the 
best two plates among the five ones obtained with the 
Schmidt telescope of the Hamburg-Bergedorf Observatory. 
The infrared comparison sequence was assembled using stars 
having a know colour index B ~V (Walker 195d) by means 
of the transformation formulae given by Allen ( 1955), taking 
care to choose only the non-reddened stars from the colour- 
colour-diagram U — B Ys B — V . 

The magnitudes of V582 Mon (KH-15D) were deter- 
mined through visual estimates with respect to its closest 
comparison stars. A stereoscopic Zeiss microscope at low and 
variable magnification was used. This method allows a com- 
parison of the measurememts obtained with different tele- 
scopes. For example, in the plates obtained with the Loiano 
reflector, the star falls in a region affected by coma, however 
the eye can perform a good estimation when comparison 
stars of comparable coma are chosen. On the basis of the 
thousands of magnitude estimations done with the above 
mentioned instruments, we estimate the mean errors in the 
magnitude estimations are: ±0.1 mag for the 0.6m Loiano 
reflector plates, ±0.05/0.1 mag for the 1.2m Asiago reflec- 
tor plates, ±0.2 and ±0.1 mag for the two 0.50/0. 40m and 
0.90/0. 65m Asiago Schmidt telescope plates respectively. 



3 DATA ANALYSIS AND RESULTS 
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Figure 2. Uppermost panel: periodogram plot of the blue light 
curve in the best sampled interval JD 2436102 (1957.725) to JD 
2438111 (1963.225). Dotted line shows the 1% false alarm signif- 
icance level, under the hypothesis of fluctuations dominated by 
Poisson statistics. The timescale corresponding to the significa- 
tive peak is 48.5 ± 0.3 days, in agreement with the phase-plots 
findings. Second panel: auto— correlations of the blue light curves, 
calculated with the DACF and ZDACF. These show a weak broad 
peak at around 48-49 days. Third panel: the first order structure 
function (SF) of the data in the same period, the SF indicates 
a drop around a timescale of 48 days. Bottom panel: the phase 
dispersion minimization (PDM) curve of the same data. The sig- 
nificant drop in the PDM gives a period of 48.4 ± 0.2 days, in 
agreement with the previous results. The results obtained with 
the different statistical methods indicate a period equal to the 
actual value. 



Part of the analysis of the pla tes taken around NGC 2264 
was published injMaffci ( 19 66al lbl). The detailed study of the 
more than 200 variable stars identified in the area of NGC 
2264 (about one square degree) is an ongoing project. At 
present, this effort has been devoted to the variable V582 
Mon (KH-15D). The infrared plates obtained with the two 
Asiago Schmidt telescopes were made available to J. A. 



Johnson, w ho analyzed the data an d published a first set 
of results in l.Iohnson fc Wim] (|2Q^. 

As shown in Tables and |5| the blue light curve is 
better sampled. However for the years 1967-70, the number 
of useful plates is nearly equal for the two colours, since after 
the positive results obtained with the first observations in 
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Figure 3. Periodogram plot of the blue light curve in the inter- 
val JD 2435394 (1955.787) to JD 2436198 (1957.988). Dotted line 
shows the 1% false alarm significance level, under the hypoth- 
esis of fluctuations dominated by Poisson statistics. The curve 
does not display any characteristic timescale or periodical fea- 
ture. Taking into account the good data sampling in 1957 (see 
Table 13 and the results displayed in Fig. |2] this implies that the 
variability/periodicity in the blue began in the first days of 1958. 
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Figure 4. Upper panel: periodogram plot of the infrared light 
curve in the interval JD 2436101 (1957.722) to JD 2438057 
(1963.077). Dotted line shows the 1% false alarm significance 
level, under the hypothesis of fluctuations dominated by Poisson 
statistics. No comparable indication of a periodicity is seen in the 
infrared light curve periodogram. Lower panel: periodogram plot 
of the infrared light curve in the period 1967-1970. In this inter- 
val, the characteristic periodicity (48.4 day) is clearly observed 
also in the infrared data. 



the infrared, PM decided to take both blue and IR images 
for every available observing night. 

In the blue, V582 Mon appears with a constant mag- 
nitude from 1955 to 1957 (Fig. 0. The first clear sign of a 
significant weakening (0.4 mag) was observed on January 8, 
1958. In the infrared, the star appeared with constant mag- 
nitude in the subsequent 4 years, and a first sure weakening 
(0.5 mag) is observed on November 16, 1961 (Fig.Ql. 



3.1 Temporal Analysis 

In addition to the classical phase analysis, the temporal 
analysis of the data was carried out using well tested and 



robust techniques suitable for an unevenly sampled time se- 
ries: the Lomb-Scargle periodogram (LSP), the first order 
structure function (SF), the discrete auto-correlation func- 
tion (DACF), and the phase dispersion minimization (PDM) 
method . These methods allow to investigate quantitatively 
the statistical temporal structure of variability, the auto- 
correlations, and the existence of characteristic time scales 
and periodicity. 

The SF is equivalent to a power spectrum calculated 
in the time domain. It provides information on the type of 
variability and the range of the characteri stic time scales 
that contr ibute to th e fluctuations (see, e.g. iRutman Ill978l : 
[Simonct ti. Cordes. fc Heeschen 1985}) . Steep drops in the SF 
indicate small variance and provide the signature of possi- 
ble characteristic time scales, while the slope in the log-log 
representation gives the value of the power law spectrum 
index, and thus the nature of the variability. The DACF 
allows to investigate the level of auto-correlation in un- 
evenly sampled data sets (se e, e.g. lEdelson fc Kroli'k1ll988l : 
iHufnagel fc Breeman Il992h without any interpolation. The 
number of points per time bin can vary greatly in the DACF, 
but data bins by equal population can be built together with 
Montecarlo estimations for peaks and uncertainties as i n 
the Fisher ^-transformed DACF (ZDACF. lAlexander lim^ . 
The LSP is a technique analogous to the Fourier analysis for 
discrete unevenly sampled data trains, useful to detect the 
strength of harmonic com ponents with a certain angular fre- 
quen cy ui = 2%f (see, e.g. IScargle ll98llHorne fc Baliunas I 
| l986h. Firiallv.t he PDM method (see, e.g. ISteUingweril978l : 
l,Turkcvich ll971^ seeks to minimize the variance of the data 
at constant phase with respect to the mean value of the 
light curve. The PDM method has no preference for a par- 
ticular periodical shape, it incorporates all the data directly 
into the test statistic and it is thus well suited to small and 
randomly spaced samples. A value is statistically significant 
when the PDM drops near zero. 

In Fig. H we plot the LSP, DACF and ZDACF, SF, and 
PDM for the blue light curve data in the best sampled inter- 
val JD 2436102 (1957.725) to JD 2438111 (1963.225). The 
timescale corresponding to the significative peak is 48.5±0.3 
days for the LSP. The DACF and ZDACF also show a weak 
(value ^ 0.3) peak for a time lag at around 48-49 days. Both 
the SF and the PDM show a drop around the same value. 

The apparently stable periodicity of the fiuctuation is 
first observed in 1958. The period is in good agreement with 
that found in recent years with well sampled data. On the 
other hand the blue light curve does not display any char- 
acteristic timescale or periodical feature before 1958 (Fig. 
EJ, reinforcing the evidence for a beginning of the eclipses 
around the first months of 1958. 

The temporal analysis of the infrared light curve for 
1957-1963 does not show any significant feature above sta- 
tistical fiuctuations in the LSP, DACF, SF, and PDM plots 
(upper panel of Fig. 2J, which indicates no corresponding 
variability and periodicity in the infrared in these years. We 
believe that the variability occurred only some years later 
(possibly about 4 years later) in the infrared based on the 
periodical behaviour in the IR displayed in the last interval 
of our light curve (lower panel of Fig.QJ, and by the B — I 
colour analysis in Par. 13. 4i 
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Figure 5. Optical blue magnitude light curve of V582 Mon 
(KH 15D) as function of the phase ip from our 1955-1970 data. 
The spectroscopic bin ary period (48.38 days) and the epoch 
2452352.26 reported bv lHerbst et alJ i20n2ft are used. All the data 
points reported in the plot are real observations, without any data 
point duplication. The eclipses appears to start at the end of the 
1950s. The periodic variability in the subsequent years up to the 
1970 seems rather sinusoidal, with slow descendent and increasing 
stages and no plateaus. 
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Figure 6. Infrared magnitude light curve of V582 Mon (KH 15D) 
as a function of the phase (j> from our 1955-1970 data. The spec- 
troscopic bi nary period (48.38 days) and the epoch 2452352.26 
reported bv lHerbst et alJ <2002l) are used. All the data points re- 
ported in the plot are real observations, without any data point 
duplication. In this plot the eclipses seems possibly to start at the 
beginning of 1960s, and the periodic variability in the subsequent 
years appear approximatively sinusoidal, with slow descendent 
and increasing stages and no evident plateaus. 



3.2 Blue Light Curve 

After the first weakening, the variability amplitude increased 
until about 1.2 mag from the end of 1967 to the beginning of 
1970. Applying the phase periodical analysis, the maximum 
probability period is 48.5 days, very similar to the ephemeris 
outlined bv lHerbst et al.l l)2nn2h . The phase light curve con- 
structed with th e observations from 1967 to 1970, adopting 
the ephemeris of iHerbst et al.l i2002l) . shows a marked reg- 
ularity. 

Figure El shows the optical blue magnitude light curve 
of V582 Mon (KH 15D) as a function of the phase for 
our data. All the data points reported in the plot are real 
observations. The eclipses appear clearly to start at the end 
of the 1950s. The periodic variability in the subsequent years 
up to the 1970 seems rather sinusoidal, with slow decreasing 
and increasing stages and no evident plateaus, contrary to 
more recent observations. A phase shift of the light curve 
minimum was found in the infrared also bvljohnson fc WinnI 




to the Ic band as mentioned above) after the first weaken- 
ing the amplitude increases, beginning in 1968, by 1.2 mag. 
The minima are coincident with those observed in the blue. 
The phase-constructed light curve shows a flat minimum in 
the years 1967-1970, similar to that observed in 1995, how- 
ever it represents a smaller fraction of the period (0.2) than 
in the more recent data (0.4). The observed maximum does 
not appear as a pronounced plateau, as seen in recent obser- 
vations. In addition, the transition between minimum and 
maximum values is much slower (0.2 against 0.07 in phase 
units); the infrared light curve is better described by a si- 
nusoidal function rather than by a step-like structure, as 
observed in 1996-2002. 

The infrared maximum appears approximatively at the 
same level as in recent years. We do not see a clear weaken- 
ing, as instead suggested recently in the infrared magnitudes 
measured in the real I Cousin band. We must however cau- 
tiously underline that the two photometric systems are not 
identical. On the other hand, the variability amplitude in- 
creased from 1 mag to 3 mag. 



3.3 Infrared Light Curve 

The same behaviour is displayed by the infrared magnitude 
light curve in Fig. |S| In the photographic infrared (similar 



3.4 Colour Indexes 

The suggested 4-year difference in starting times between 
the blue and the IR periodicity is supported by the analysis 
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Figure 7. Upper panel: the B — I colour index in function of 
the time. The plot is clearly in agreement with an infrared colour 
excess in the period 1955-1963, implying that the beginning of 
the periodical eclipsing variability in the blue was not followed 
by the infrared flux that remains rather constant in those years. 
Lower panel: the scatter plot between the B — I colour index 
and the B magnitude brightness shows the same behaviour. The 
linear regression trend lines are reported. The slope is steeper for 
the 1955-1963 data (blue/darker square symbols) with respect to 
the slope of the 1967-1970 data (azure/lighter diamond symbols) 
implying that the star becomes redder when fainter during the 
1955-1963 epoch. 

of the B — I colour indexes. For the calculation of the the 
B — I index, we selected only the data from the two bands 
which were coincident within an interval of 1 day. Figure |7| 
shows the behaviour of the colour index in time as well as 
the correlation with the blue brightness, which fact confirms 
and underlines our finding that the beginning of the infrared 
variability occurred some years later (possibly 4-years) after 
that in the blue. In fact, the B—I index shows a colour excess 
and a variability for 1955-1963, while it tends to remain 
constant in subsequent years (Fig.Qupper panel). Moreover, 
the scatter plot B — I versus B shows a steeper trend for 
the 1955-1963 data, implying that the star became redder 
while fainter during these years. We point out that there 
are several parallel B and I observations within one day, to 
support this hypotheses. Only in the subsequent years (1967- 
1970) does the behaviour of the colour index not depend on 
the source brightness, implying that the star shows the same 
periodical behaviour in the blue and infrared bands. 



bands between October 13, 1955 and March 1, 1970. This 
set of homogeneous data in two colours is unique to our 
knowledge for this period and reveals the start of the eclipse 
phase of V582 Mon in the two colours. 

The analysis of our data provides the following impor- 
tant results: 

(i) in the blue, the star showed the first magnitude vari- 
ations at the beginning of 1958; 

(ii) in the near infrared, the first magnitude variation was 
observed roughly 4 years later; 

(iii) over the above 4 year period the star showed an in- 
frared excess, at least around minima. 

(iv) the period did not change; 

(v) the shape of the light curve appears difi'erent from 
the present one, i.e. a sinusoidal functi on rather than a 
curve with flat minima and maxima, as in lKearns fc Herbsti 
||1998^ . 

A possible model to explain these characteristics and their 
evolution would be the following. Until 1958, V582 Mon 
(KH 15D) was surrounded by an edge-on disk or gaseous 
torus, which did not cause variations in the brightness. Sub- 
sequently, thin dust, able to absorb the blue radiation (but 
not the IR radiation) began to form, distri buted as a sort 
of "ba nana" as in the model proposed by iBarge fc VitonI 
i2003h . but with nearby, or extended thinner edges as indi- 
cated by the absence of a fiat maximum. Afterwards the dust 
aggregated into larger particles, and the selective absorp- 
tion became total absorption. Due to the increasing num- 
ber and dimension of the dust grains, the observed minima 
deepen. At the same time, the "banana" dust cloud (or its 
denser part) contracted, creating wider, fiattened minima 
and steeper variations of the brightnes s to reach the ac - 
tual shape and values of the light curve JWinn et alJl2004ri . 
This suggests the presence of a body or agglomerate of dust 
and/or stones less extended but denser and with sharper 
edges. 

These hypotheses imply that in the last 40 years we have 
observed the birth of a giant planet near the star, as that dis- 
covered (with other techniques) in other nearby stars (see, 
e.g. I^utler 2003; Marcy & Butler 1998). More likely, we ad- 
vance the notion of the formation of a blob of material, from 
which one or more planets could form, or that could dissolve 
or evolve in an unforeesable way. This is supported by the 
fact that V582 Mon (KH 15D) is a pre-main sequence star, 
i.e. in formation. Our data does not exclude other models. 
In any case, any proposed model must take into account the 
non- variability prior to 1958, the progressive increase of the 
absorption, the interval of 4 years separating the beginning 
of the blue and infrared variability and above all, the trans- 
formation from a sinusoidal light curve to an actual step-like 
structure characterized by fiat maxima and minima. We be- 
lieve that the evolution of this unique and peculiar stellar 
object must be followed in the future with several different 
instruments, and in particular with the same photometric 
and spectroscopic techniques used in the past. 



4 DISCUSSION AND CONCLUSIONS 

We have presented the results of a long-term series of ob- 
servations of the variable star V582 Mon (KH 15D) in the 
cluster NGC 2264, performed in the blue and near-infrared 
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